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Abstract 

We present a study for the pion decay constant f n in the quenched approxi- 
mation to lattice QCD with the Kogut-Susskind (KS) quark action, with the 
emphasis given to the renormalization problems. Numerical simulations are 
carried out at the couplings = 6.0 and 6.2 on 32 3 x 64 and 48 3 x 64 lat- 
tices, respectively. The pion decay constant is evaluated for all KS flavors via 
gauge invariant and non-invariant axial vector currents with the renormaliza- 
tion constants calculated by both non-perturbative method and perturbation 
theory. We obtain f n = 89(6) MeV in the continuum limit as the best value 
using the partially conserved axial vector current, which requires no renor- 
malization. From a study for the other KS flavors we find that the results 
obtained with the non-perturbative renormalization constants are well con- 
vergent among the KS flavors in the continuum limit, confirming restoration 
of SU(4)a flavor symmetry, while perturbative renormalization still leaves an 
apparent flavor breaking effect even in the continuum limit. 

PACS number(s): 12.38.Gc, 12.39.Hg, 13.20.He, 14.40.Nd 
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I. INTRODUCTION 



In recent large-scale simulations of lattice QCD, statistical errors of physical quantities 
have become quite small. Indeed, for some hadronic matrix elements the precision has 
been so high that we cannot ignore uncertainties coming from the renormalization factor of 
lattice operators. Thus it has become increasingly important to reduce uncertainties from 
this source. 

Renormalization factors can be evaluated in perturbation theory. Pushing the calculation 
beyond the one- loop level is usually difficult, however, hence uncertainties arising from 
higher-order corrections remain. We expect this problem of higher-order uncertainties to 
profit fully from a non-perturbative treatment. A non-perturbative method for calculating 
renormalization factors was proposed JIJ , and has been applied to the quark mass || , decay 
constants || and four-fermion operators 0, with the Wilson and the clover quark actions 
and to the quark mass with the Kogut-Susskind quark action ||. 

An important point to check with non-perturbatively calculated renormalization fac- 
tors is their reliability and the degree of improvement achieved in the final physical results. 
For this purpose the pion decay constant is perhaps the best choice because the reference 
experimental value is known to a high precision. A verification that non-perturbative deter- 
mination works for simple quark bilinear operators is a first step to ensure validity of more 
general applications to four-quark or other operators. 

In this work, the pion decay constant is examined with the Kogut-Susskind (KS) quark 
action via gauge invariant and non-invariant operators using all KS flavors. The KS action 
has the well-known feature that SU(4)a flavor symmetry is broken down to U(1)a subgroup 
at a finite lattice spacing. We orient our study mainly toward the following two points 
provided by this feature. First, due to the remaining U(1)a symmetry, the renormalization 
constant for the corresponding axial vector current equals exactly unity, and hence the pion 
decay constant calculated in this channel receives no renormalization. This makes it possible 
to attain a high-precision calculation of the pion decay constant without uncertainties from 
renormalization. Second, we can calculate the pion decay constant using axial vector currents 
in the other KS flavor channels. Symmetry is broken in the decay constants at a finite lattice 
spacing, but restoration is expected in the continuum limit. Such restoration of full flavor 
symmetry has been previously examined for pion mass (Bl[7|. Here we extend the study to the 
pion decay constant, the new feature being the necessity of renormalization constants. This 
can be used to investigate the reliability of non-perturbative methods for the calculation of 
renormalization factors, compared to perturbative treatments. We also compare the results 
obtained with gauge invariant operators to those with non-invariant ones. 

The paper is organized as follows. In Sec. [TJ we establish our notations and formalism. 
The method employed for our calculations is explained in Sec. |T|, followed by discussion of 
perturbative and non-perturbative renormalization factors in Sec. |TV| We summarize the 
simulation details in Sec. [V], and present the results on the chiral and continuum extrapola- 
tions in Sees. [VI| and |V1J] , We close with a brief conclusion in Sec. |V111| . 
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II. FORMULATIONS 



A. Kogut-Susskind quark action 

The Kogut-Susskind (KS) quark action is defined in terms of one-component fermion 
fields x( n ) an d x( n ) on a lattice whose site is labelled by = 0, 1, 2, . . . , L — 1, 



S, 



-X(n + p,)Ul(n)x(n)) + m q x(n)x(n) 



(1) 



where m q is the bare quark mass and t]^ = (— is the KS sign factor. Color sums 
are assumed for simplicity. Dividing the lattice into 2 4 -hypercubes which are labelled by 
Xfj, — 0, 2, 4, . . . , L — 2, and whose corners are specified by a four- vector A with — or 1, 
we introduce sixteen-component fields 



1_ 1 

a( x ) = iX{x+ A ), 4>a{x) = -x(x + A). 



In terms of these hypercubic fields, the action (HD is rewritten as 



(2) 



Sf = (2a) 4 £ 



x,AB 



+a0 A (a;)(75 ® i^)AB^U AB (j) B (x)) 
+m q 4> A {x){I ® I)abU A b{x,x)(/) B (x) 



(3) 



Here a hypercube matrix referring to the Dirac spinor 75 = 7^ • • ^f 4 and the KS flavor 
£f = 7i Fl • • • 74 Fi is defined by 



(75 ® £f)ab = ^Tr (7a757b7f) 



(4) 



and the lattice derivatives are given by 



V^UAB^Bix) = — (U AB (x,y)(f)B(y)\y=x+2p, 
-U AB {x,y)(j)B(y)\y=x-2p) 



&pUab4>b{x) = ^(UAB(x,y)4> B (y)\y=x+2fL 

+UAB(x,y)<f>B(y)\v=x-2fi 

-2U AB (x,x)(j> B (y)), 



(5) 



(6) 



where U A B(x,y) is the average of ordered products of gauge link variables over the shortest 
paths from x + A to y + B. 
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A gauge-invariant meson operator with zero spatial momentum is defined in this hyper- 
cubic notation by 

°s( x d =J2H^a( x )(1s ® £f)abUab(x,x)(I>b(x). (7) 

x AB 

For instance, O5 = A^, ir, and pk respectively for 75 = 7^75, 75, and 7^. Here we have 
2 4 x 2 4 = 256 operators, which are classified into irreducible representations || in terms of 
dn = Sfj,- (mod2). 

The form of the action (|3]) shows that the flavor-mixing term (75 £g> £^£5) breaks SU(4)a 
flavor symmetry down to U(1)a subgroup for the flavor channel £5 at a finite lattice spacing. 
A lattice analog of the PCAC relation holds in the £5 channel corresponding to U(1)a 
symmetry: 

V M 4(x) = 2m,7r 5 (x), (8) 

where the superscript 5 refers to £5. On the other hand, there appear additional terms in 
the PCAC relation for other channels, which vanish only in the continuum limit. 

B. Pion decay constant 

The pion decay constant is defined in the continuum theory by 

V2Um w = (0\u 7a5 d\n + (p = 0)>. (9) 

We adopt the normalization f^ exp ^ — 93 MeV. If we use the PCAC relation, this may be 
rewritten as 

V2f w ml = (m u + m d )(0\u l5 d\7r + (p = 0)). (10) 

The lattice pion decay constant for the KS flavor £p is defined by 

V2f^ = (0\A^ F (p = 0)), (11) 

In the £5 channel where the PCAC relation (§) holds, we may use an alternative formula 
corresponding to Eq. (|T0|) : 

v^/f^K) 2 = 2m g (0|7rV(p = 0)), (12) 

where we have added the superscript (P) to distinguish explicitly the pion decay constant 
obtained with a pion operator from that with an axial vector current. 

III. EXTRACTION OF PION DECAY CONSTANT 

We employ the wall source technique to enhance signals || . The meson operator for the 
wall source at the origin is defined by 
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(13) 



where we assume that gauge configurations are fixed to some gauge. The matrix elements 
appearing in the definition of the pion decay constant are extracted from the large-time 
behavior of the correlation function at zero spatial momentum: 

<o£(*)*£(o)> 

~ <^> t Hexp(-m^) ± exp(— mf (T - t))], 
f + sign for O s = it, tc w , 

\ - sign for O s = At, { > 

where is pion mass common to the three cases. Here we extend the time slice of meson 
operator defined at x 4 = 0, 2, 4, . . . , T — 2 to have t — 0, 1, 2, . . . , T — 1 extensions with the 
temporal lattice size T. Note that there is no mixing between a t — ±1 states in this case 
with time-extended meson operators. The amplitude Cq s ^ w can be written up to an overall 
sign factor as 

rF (0|Of|^(p = 0))( 7 r^=0)|^|0) 

with V, the spatial lattice volume. Using the amplitude of the correlation functions with 
the axial vector current (Os = At) and the pion operator for the wall source (Os = nw), 
the pion decay constant is calculated as 

fmP fC* ' 

where the pion mass obtained by the correlation function with the pion operator (Os = 7r) 
is used in this work. For comparison, the gauge non-invariant axial vector current and pion 
operator to obtain the amplitude and pion mass, respectively, 

A'/(x±) = EE^)(747s ® Zf)abMx), (17) 

x AB 

k >F (xi)=J2Y.Mx)(15®Zf)abMx), (18) 

x AB 

is also examined. Alternatively, an extraction of the decay constant from the pion operator 
(Of = 7r 5 ) requires the combination given by 



f (P)5 _ V 2 ^* °™w (19) 



TTW^W 
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IV. RENORMALIZATION 



A. General considerations 

Renormalization is necessary to extract the physical pion decay constant from the lattice 
calculations. This procedure is made for each flavor in the case of the KS action. It is 
expected that the renormalization eliminates the KS flavor dependence in a way that the 
decay constant calculated for various KS flavors takes a unique value in the continuum limit. 

Let us define a multiplicative renormalization constant for the lattice axial vector 
current A^| lat through 



A 



phys 



lat 



(20) 



According to the definition (|TTD the pion decay constant calculated with the axial vector 
current is renormalized as 



J TT 



phys 



r/F fF 
A J TT 



lat 



(21) 



As a special case, we have 



L A 



(22) 



in the £5 channel due to the lattice PCAC relation (El). Thus the pion decay constant can 
be calculated with out any uncertainties of renormalization in this channel, while the other 
channels can be used to check the reliability of renormalization constants by examining 
the expected convergence of the renormalized pion decay constants to a single value in the 
continuum limit. 

The decay constant defined with the pion operator (pT2|) is renormalized as 



J TT 



(P)5| 



I phys 



(P)5| 



I lat' 



(23) 



where Z m is the renormalization constant for quark mass. Using the identities Z m — \jZ\ 
and Zg = Zp, where the superscript / refers to the KS flavor for a unit matrix, we find that 
this relation is identical to 



■J TT 



(P)5| 



I phys 



■J TT 



(P)5| 



I lat' 



(24) 



which is equivalent to Eq. 



B. Perturbative and non-perturbative renormalization factors for axial vector 

currents 

We employ two sets of the renormalization factor Z\ for the KS axial vector current. One 
of them is perturbatively calculated at one-loop order |T0|. We apply tadpole improvement to 
the axial vector current operator using the fourth root of plaquette as the tadpole factor, and 
evaluate the renormalization constants with the tadpole-improved MS coupling at q* — 1/a. 
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The other is non-perturbatively evaluated with the regularization independent (RI) scheme 
of Ref. IH, which was developed for the Wilson and clover actions. In the RI scheme, the 
renormalization factor is obtained from the amputated Green function in momentum space 

r£» = 5(p)- 1 (O|0(p)Of0(p)|O)5(p)- 1 , (25) 

where the quark two-point function is defined by S(p) = (0\<f)(p)<f)(p)\0), and the momentum 
of the hypercubic field <p(p) takes values of the form p^ = 27rn^/(aL) with — L/4 < < 
L/4 — 1. The renormalization condition imposed upon Tq s (p) is given by 

Z ( Os )F (p)Mp) = mK^Osip)], (26) 

where (j^Osf = (7s ® £f) 1S the projector onto the tree-level amputated Green function. 
The wave function renormalization constant is calculated by imposing the condition 
Z v (p) = 1 for the conserved vector current for (7^ ® /). The relation between the overall 
renormalization constant Z^ appearing in Eq. (^) and Z^ is simply 

Z F = l/zf )F , (27) 

because the continuum axial vector current is not renormalized. 

The calculations for the non-perturbative renormalization constants were carried out in 
quenched QCD in our previous publication ||. The results for the scalar and pseudoscalar 
operators have been used in our analysis of light quark masses for the KS quark action in 
quenched QCD ||. Here we use them for the axial vector renormalization factors. 

The calculational parameters are summarized in Table |. We evaluate the Green func- 
tion ( p5|) for 15 momenta in the range 0.038533 < (pa) 2 < 1.9277 using quark propagators 



evaluated with a source in a momentum eigenstate. In Fig. [I] we present the renormalization 
constant for both vector and axial vector currents, respectively denoted by Z v R '^ F (p) and 

%A (p)i i n the chiral limit. 

A practically important issue with the non-perturbative method employed here is the 
choice of the momentum at which the renormalization factors are evaluated. In general the 
momentum should satisfy Aqcd C |) < 0(a _1 ) in order to keep under control the non- 
perturbative hadronization effects and the discretization error on the lattice. Since these 
effects appear as p-dependences of renormalization factors, we should avoid the range where 
a momentum dependence is visible. Another point to consider is the relation Zy F5 (p) = 
Z^ F (p) with the superscript F5 referring £^£5, which we would expect to hold for all 
momenta p in the chiral limit due to U(1)a chiral symmetry of the KS quark action. 

For (3 = 6.2 Fig. [I] shows that these two requirements are satisfied for p 2 > 5 GeV 2 , 
which corresponds to (pa) 2 > 0.5. In order to satisfy p 0(a^ 1 ), we take (pa) 2 = 1.0024 
(p 2 = 7.0392 GeV 2 in physical units) to calculate the renormalization factors used for the 
pion decay constant. The same value of lattice momentum (pa) 2 = 1.0024 is chosen for (3 = 
6.0, which corresponds to p 2 = 3.5428 GeV 2 . The numerical values of the renormalization 
factors are summarized in Table 0. 
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V. DETAILS OF SIMULATION 



A. Simulation parameters 

We carry out our calculations in quenched QCD using the standard plaquette action 
for gluons. As we summarize in Table [TTT| , numerical simulations are carried out at (3 = 
Q/g 2 = 6.0 and 6.2 on 32 3 x 64 and 48 3 x 64 lattices, respectively. Gauge configurations are 
generated with the five-hit pseudo-heatbath algorithm, and hadron correlation functions are 
calculated on 100(60) configurations separated by 2000 sweeps at (3 = 6.0(6.2). 

Gauge configurations are fixed to the Landau gauge through maximization of 

F L = 2>(tf M (n) + E>J(n)). (28) 

This is realized by iterating the steepest descent method for the first 2000 steps and the 
over-relaxation method for the subsequent 3000 steps until the condition 

A = ^7Etr(G[HG L (n))<10- 14 (29) 

is satisfied, where V is the lattice volume and 

Gl = Y1 \( u ^( n ) ~ U fi(n - ft) - h.c. - trace). (30) 

We take three values for quark mass, m q a = 0.030, 0.020, 0.010 at (3 = 6.0 and 0.023, 
0.015, 0.008 at (3 = 6.2. Quark propagators are evaluated for 16 types of wall sources, each 
corresponding to a corner of a hypercube, defined by 

J2 d ab(x,v) J2 G bc(v,z) = J2 5 *z S ac, (31) 

y,B z z 

where Dab(x, y) is the quark matrix for the KS action. We solve the equation independently 
for each C by the conjugate gradient method with the stopping condition 

||remnant vector|| 2 < 10~ 5 . (32) 

The 16 quark propagators are combined to construct the 16 meson correlation functions in 
the KS flavor basis specified by the hypercube matrix £p. Averages are taken of the meson 
correlation functions over 2 3 ways of choosing the spatial origin of hypercubes on the lattice. 
We also average them over all states belonging to the same irreducible representation |J. 

B. Fitting procedure 

In fitting the meson correlation function C(t) to the asymptotic form C fit (t) for an 
extraction of the mass and amplitude, we symmetrize the correlator at t and T — t, and 
carry out a standard correlated fit minimizing 
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X 2 = J2AC(t)X-\t,t')AC(t'), 
t,v 



(33) 



where 



Z(t,t') = {C(t)C(t'))-{C(t)){C(l/)) 



(34) 



is the covariance matrix of the correlator and AC(t) = C(t) — C fit (t). The fitting range 
t = t min , . . . ,t max is chosen by fixing t max = T/2 and varying t min so that x 2 /iV DF takes a 
value near unity, where A^f is the degree of freedom of the fit. Finally, errors in this work 
are estimated by the single elimination jackknife procedure. 



We check the validity of the asymptotic form of the mesonic correlation function (|14"|) 
which is based on the assumption of a single pole dominance by an inspection of the effective 
mass. Typical results for the effective mass extracted from the correlators (n F (t)n^(0)) and 
(A F (t)vr^(O)) are compared in Fig. ||. We observe a wide plateau and an expected agreement 
of the effective masses from the two correlation functions. We then find no problem in fitting 
these correlation functions by a single pole. 

The situation is different for the wall-to-wall correlation function (Tr^(^)vr^(O)), particu- 
larly at P = 6.2. As we show in Fig. |3|, the effective mass for {^wi^) 7l w(^)) does no ^ reacn a 
plateau at (3 = 6.2 even at t ~ T/2, and agreement with the effective mass of (7r F (t)7r^(0)) 
is not seen. This behavior is most likely caused by a lack of sufficient temporal size of the 
lattice, and poses a practical problem of how one extracts the wall-to-wall amplitude C F w7Tw 
which is needed in Eq. (|16]) to calculate the pion decay constant. 

To solve this problem, we perform a double pole fit for (^wi^) 7l w(^)) gi ven by 



Ideally one likes to make a fit with four parameters C F w7Tw , m F , C F g, and m F ^. This fit, 
however, is quite unstable because the fitting function consists of a sum of two exponentials 
with not much different masses m F and m^. Therefore, we fix the pion mass parameter m F 
to that obtained from (7r F (t)7r^(0)). 

As we now can no longer compare the effective pion mass for (Tr^(t)vr^(O)) to that for 
(n F (t)n F (0)) , we present a typical comparison of the amplitudes, extracted with the fitting 
range from t to T/2 with the single and double pole fits, in Fig. ^. We also compare xV^df 
for the two fits in Fig. [5]. From these figures, we consider that the double pole fit provides 
a good determination of the amplitude C F of the pion to the wall operator with a wide 
plateau of the amplitude and a reasonable value of xV^df ~ 0{1). 

A possible interpretation for the dominant source of contamination to the wall-to-wall 
correlation function is an unbound quark-antiquark pair. Such a state can contribute since 
gauge configurations are fixed to the Landau gauge. In Fig. ^|we plot the value of the second 



C. Wall-to-wall amplitude 



<7r£(t)7r£(0)> 
~ ^™(^[exp(-m: t) + exp(-m: (T - t))] 
+ Cgg [exp(-m£t) + exp(-m^(T - t))}. 



(35) 
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pole mass m qq as a function of quark mass. The fact that the results depend little on the 
KS flavor of the meson operator is consistent with this interpretation. In the chiral limit 
one obtains m qq ~ 2 x 440 MeV, which is a reasonable value for a constituent quark mass. 

Finally, we summarize the fitting ranges £ m j n common for all flavors and xV^df for our 
global fits in Table [IV]. Here, we have used the alternative fitting range of the wall-to-wall 
correlation function to improve the fitting quality for £p = £4, because the common fitting 
range does not give a satisfactory result [TT| caused by worse fitting. 



VI. CHIRAL BEHAVIOR 
A. Pion masses 

We show values of (m^a) 2 as a function of m q a in Fig. [7|. Pions for the 16 KS flavors 
are classified into 8 irreducible representations. These consist of four 1-dimensional rep- 
resentations given by £ 5 , £ 4 £ 5 , £4, I and four 3-dimensional representations given by ^£5, 
6c£4) ik^li (k,£ = 1,2,3; k < £). We observe very clearly in Fig. [7] that these irreducible 
representations form a degeneracy pattern specified by 

6, (£4, £fc), /• (36) 

This pattern was observed long time ago in Ref . B . A theoretical explanation based on the 



effective chiral Lagrangian analysis for KS quark action was provided recently in Ref. [[L^ . 

Another notable feature in Fig. ^ is a linear behavior of pion masses as a function of 
quark mass from the correlation function with the gauge invariant pion operator. With a 
linear extrapolation we observe a non-vanishing value at m q a = in channels other than £5 
for which U(1)a symmetry holds. The gauge non-invariant case, not presented in the figure 
but in Table [V] for the numerical values, also shows almost the same result as in Fig. [7[ 

The chiral behavior of p meson mass for various KS flavors is shown in Fig. |^. We find 
the difference of masses among various flavor channels to be small, less than 1% even in the 
chiral limit obtained by a linear extrapolation. We therefore choose the p meson mass in 
the flavor channel (7^ <8> £&), for which the p meson operator is local, to set the scale using 
the experimental value m p = 770 MeV. We then find that a -1 = 1.92(2) GeV for (3 = 6.0 
and a- 1 = 2.70(5) GeV for (3 = 6.2. 



B. Pion decay constant 

In Fig. [J] we illustrate the chiral behavior of the bare pion decay constants calculated 
with Eq. ([TB]). As with the case for pion masses, we use a linear extrapolation toward the 
chiral limit. 

The pion decay constants obtained for eight irreducible representations again form a 
degeneracy pattern, which, however, is different from that for pion masses. This is due to 
the fact that the pattern for the decay constant reflects the distance of the axial vector 
current operator rather than that of the pion operator: the two operators differ because of 
the the Dirac factor, 7475 for the axial vector current and 75 for the pion. We also observe 
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that the KS flavor dependence of the decay constant is much larger for the gauge invariant 
operators than that for the non-invariant ones. In contrast to the case of mass, for which no 
renormalization is required and lattice symmetry group controls, the pattern for pion decay 
constants mainly comes from the insertion of gauge link variables, which is roughly written 
as relation between the continuum and lattice axial vector currents: 

1 ,\ 4 , F 



A»L nt ~[^Ua)) AXv (37) 

Here d is the distance of the axial vector current operator for the gauge invariant case, while 
the non-invariant operator corresponds to d — 0. 



We show the decay constants after renormalization in Figs. (TTJ and IT]. With the use of 



perturbative renormalization constants (Fig. ITUp, the discrepancy among different KS flavor 



channels becomes smaller toward the continuum. The reduction of the discrepancy, however, 
is significantly more dramatic with the use of non-perturbative renormalization constants as 
shown in Fig. |ll[ In particular, the large difference among bare results obtained with gauge 
invariant operators almost disappears. 

The numerical values for pion decay constants are collected in Tables |V"1|^V111| . In contrast 
to the case of pion mass, there is no flavor channel to give the same results for the gauge 
invariant and non-invariant case, because the simultaneous local channel does not exist for 
the axial vector current and the pion operator both appearing in the calculation of the pion 
decay constant. 



VII. CONTINUUM EXTRAPOLATION 

In Fig. UJ], we present a-dependence of (mf ) 2 quadratically extrapolated to m q a = 0, 
according to 0(a 2 ) scale violation expected for the KS quark action. We observe clear 
evidence that the non-zero values of (m^) 2 for the non-Nambu-Goldstone channels vanish 
as a 2 toward the continuum limit, supporting the restoration of full flavor symmetry of the 
KS action. 

The continuum extrapolation of the pion decay constant, renormalized perturbatively or 
non-perturbatively, is shown in Fig. [13] as a function of a 2 . In this figure with an enlarged 



vertical scale as compared to Figs. |IJJ and |ll], we observe a general trend that the difference of 
values among various KS flavors becomes smaller toward the continuum limit. In particular, 
for non-perturbatively renormalized decay constants the central values in the continuum 
limit agree within a 2% accuracy, which is well below the statistical errors of 5-10%. On the 
other hand, the convergence is worse for the perturbatively renormalized decay constants. 
The spread in the continuum limit is 3-4%, which is roughly the magnitude of uncertainty 
one expects from higher-order corrections in the renormalization factors. We consider that 
these results provide evidence for both restoration of SU(4)a flavor symmetry of the KS 
action in the continuum limit and the effectiveness of the non-perturbatively evaluated 
renormalization constants. 



The values of pion mass squared for various KS flavors are listed in Table [TX|, and those 



for pion decay constants are collected in Tables |X| and X|. As our best value for the decay 



constant, we take f n = 89(6) MeV obtained with the gauge invariant axial vector current 
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in the £5 channel which requires no renormalization. This value is compared with the 



experiment 92.4(3) MeV [13|. Possible quenching errors are not visible within the statistical 
error of 6 MeV. 

Let us recall that the decay constant in the £5 channel can also be calculated from the 
pion operator using Eqs. fll9"|) and ([?|). Results are added in the bottom lines of Table [X| 
(and |Xj for convenience of readers), which show reasonable agreement with those from the 
axial vector current in the £5 channel, as expected. 



VIII. CONCLUSION 

In this article we have presented an analysis of the pion decay constant in quenched 
QCD using the Kogut-Susskind quark action. Our best estimate for the decay constant in 
the continuum limit is 89(6) MeV, which is obtained with the gauge invariant axial vector 
current which respects U(1)a symmetry. 

We have carried out a detailed comparison of perturbative and non-perturbative axial 
vector renormalization treatments. We conclude that the non-perturbative renormalization 
factors efficiently eliminate the flavor breaking effect in the decay constant in the contin- 
uum limit, while an apparent flavor- dependent difference still remains with the perturbative 
factors. 
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FIG. 1. Non-perturbative renormalization constants for vector Zy(p) and axial vector currents 
Z^{p) examined with gauge invariant current for (a) (3 = 6.0 and (b) (3 = 6.2, and by gauge 
non-invariant current for (c) (3 = 6.0 and (d) j3 = 6.2. 
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FIG. 2. Typical comparison of the effective mass of axial vector current (open symbols) and 
pion (filled symbols) correlation functions with its global mass (horizontal lines) for £p = £5 at (a) 
= 6.0 and (b) = 6.2. Circles, squares and diamonds refer to the quark masses in descending 
order at each coupling respectively. Note that result does not depend on gauge invariance of the 
operator in the case using local operator like this case. 
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FIG. 3. Typical comparison of the effective masses of the wall-to-wall correlation function (open 
symbols) and that for pion correlation function (filled symbols) with its global mass (horizontal 
lines) for £p = £5 at (a) = 6.0 and (b) = 6.2. 
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FIG. 4. Typical comparison of of the amplitude of wall-to-wall pion correlation function for 
£f = £5 at (3 = 6.2 obtained by (a) the single pole fit and (b) the double pole fit. 




FIG. 5. Typical comparison of x 2 /Ndf for (a) the single pole fitting and (b) the double pole 
fitting of wall-to-wall pion correlation function for £p = £5 at (3 = 6.2. 
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FIG. 6. Chiral behavior of the alternative pole masses appearing in the wall-to-wall correlation 
function at (a) (3 = 6.0 and (b) (5 = 6.2. Shape of symbols refer to the distance of the opera- 
tor. Some symbols (square, diamond and up-triangle) denotes two flavors; the former one refers 
time- local operators (filled symbols) including flavor £5, and the latter one refers time-separated 
operators (open symbols). 
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FIG. 7. Chiral behavior of pion masses obtained with the gauge invariant pion operators at (a) 
j3 = 6.0 and (b) = 6.2. Shape of symbols refer to the distance of the operator. Some symbols 
denotes two flavors; the former one refers time-local operators (filled symbols) including flavor £5, 
and the latter one refers time-separated operators (open symbols). For gauge non-invariant result, 
see Table 0. 
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FIG. 8. Chiral behavior of gauge-invariant p meson masses at (a) [3 = 6.0 and (b) (3 = 6.2. 
Symbols refer to the distance of the operator, Some symbols denotes four flavors; the first two 
flavors refer time- local operators (filled symbols), and the last two flavors refer time-separated 
operators (open symbols). 
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FIG. 9. Chiral behavior of the bare pion decay constants obtained by gauge invariant axial 
vector current for (a) /3 = 6.0 and (b) (3 = 6.2, and by gauge non-invariant current for (c) (3 = 6.0 
and (d) (5 = 6.2. Omitted legends in the top two figures are same as that in the bottom figures. 
Shape of symbols refer to the distance of the operator. Some symbols denotes two flavors; the 
former one refers time-separated operators (filled symbols) including flavor £5, and the latter one 
refers time- local operators (open symbols). 
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FIG. 10. Pion decay constant renormalized by one- loop perturbative renormalization factor 
Z^' F f£ . Figures (a)-(d) correspond to those in Fig. U. 
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FIG. 11. Pion decay constant renormalized by non-perturbative renormalization factor 
Z^ F . Figures (a)-(d) correspond to those in Fig. U. 
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FIG. 12. Continuum limit of pion mass squared. Symbols are same as those in Fig. f?|. 
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FIG. 13. Continuum limit of renormalized pion decay constants. Results obtained with per- 
turbative renormalization factors for (a) gauge-invariant and (b) non-invariant operator, and those 
with non-perturbative factors for (c) gauge-invariant and (d) non-invariant operator are shown. 
Symbols are same as those in Figs. |9|— |TT|. 
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TABLES 



TABLE I. Calculation parameters 


for evaluation of non 


-perturbative renormalization constants. 


P 


L 3 x T 


m q a 


a-^GeV) 


#conf. 


6.0 


32 3 x 32 


0.010, 0.020, 0.030 


1.88(4) 


30 


6.2 


32 3 x 32 


0.008, 0.015, 0.023 


2.65(9) 


30 


TABLE II. 


Renormalization constants used for 


renormalizing pion decay constants. 


(a) = 6.0 




Perturbative 


Non-pert urbative 


Operator 


Gauge inv. 


Non-inv. 


Gauge inv. 


Non-inv. 


(7475 ® £5) 


1 


0.8917 


1 


0.85019(7) 


(7475 O £fc£s) 


1.1436 


0.8547 


1.2008(1) 


0.8527(1) 


(7475 ® £fc£l) 


1.3749 


0.8556 


1.4799(2) 


0.8656(1) 


(7475 ® U) 


1.4950 


0.8569 


1.8242(3) 


0.8736(2) 


(7475 ® £46) 


0.7908 


0.7908 


0.7976(2) 


0.7976(2) 


(7475 ® &£m) 


0.9294 


0.8277 


0.9860(1) 


0.8508(1) 


(7475 ® 6c) 


1.1440 


0.8550 


1.2294(3) 


0.8767(2) 


(7475 ® I) 


1.3837 


0.8605 


1.5145(5) 


0.8835(2) 




(b) /3 = 6.2 




Perturbative 


Non-pert urbative 


Operator 


Gauge inv. 


Non-inv. 


Gauge inv. 


Non-inv. 


(7475 ® 6) 


1 


0.8917 


1 


0.86430(7) 


(7475 O £fc&) 


1.1338 


0.8643 


1.1783(1) 


0.86363(7) 


(7475 O £fc£i) 


1.3434 


0.8651 


1.4221(2) 


0.8739(2) 


(7475 ® £4) 


1.4567 


0.8663 


1.7164(3) 


0.8803(2) 


(7475 O 


0.8065 


0.8065 


0.8136(1) 


0.8136(1) 


(7475 ® &£m) 


0.9369 


0.8401 


0.9838(1) 


0.8600(1) 


(7475 ® £fc) 


1.1342 


0.8646 


1.1999(1) 


0.8825(2) 


(7475 ® I) 


1.3508 


0.8696 


1.4472(4) 


0.8882(2) 



TABLE III. Calculation parameters of our simulation. 



p 


L 3 x T 


m q a 


a-^GeV) 


#conf. 


6.0 


32 3 x 64 


0.010, 0.020, 0.030 


1.92(2) 


100 


6.2 


48 3 x 64 


0.008, 0.015, 0.023 


2.70(5) 


60 



23 



TABLE IV. Minimum time slice t m i n common for all flavors except for £p = £4 in the paren- 
thesis (See text for reason), and xV-^df of global fits for the local channel. 



(Af(t)TT^(O)) (vr^(t)Tr^(O)) (vr^(t)vr^(O)) <pf (t)p^(O)) 

(3 m q a t mill x 2 /Ndf t min x 2 /Ndf t mia X 2 /Ndf t min X 2 /Npf 

6.0 0.030 17 1.37 17 1.24 14 1.07 18 1.34 

0.020 17 1.05 17 0.95 15 1.27 17 0.87 

0.010 16 0.97 16 0.99 15 0.83 15 1.27 

6.2 0.023 17 1.40 24 1.06 17(16) 0.92 23 0.79 

0.015 16 1.07 23 0.85 19(19) 0.97 23 0.94 

0.008 15 1.33 22 0.99 19(20) 1.21 22 0.58 



TABLE V. Pion mass squared (m^) 2 in lattice units. Note that the correlation function with 
the local pion operator in the £5 channel gives exactly the same results for the gauge invariant and 
non-invariant case. 

(a) (3 = 6.0 

Gauge invariant Non-invariant 



Operator m q a = 0.030 m q a = 0.020 m q a = 0.010 m q a — * m q a = 0.030 m q a = 0.020 m q a = 0.010 m q a — > 



(75 t 


5&) 


0.1687(3) 


0.1129(3) 


0.0575(2) 


0.0018(2) 










(75 i 




0.2077(4) 


0.1454(4) 


0.0846(4) 


0.0228(4) 


0.2077(4) 


0.1454(4) 


0.0846(4) 


0.0228(4) 


(75 i 




0.2194(4) 


0.1561(5) 


0.0946(6) 


0.0317(6) 


0.2194(4) 


0.1562(5) 


0.0947(6) 


0.0317(5) 


(75 1 


3 €4) 


0.2260(5) 


0.1630(6) 


0.1023(8) 


0.0396(9) 


0.2261(5) 


0.1630(6) 


0.1024(8) 


0.0396(9) 


(75 i 




0.2086(5) 


0.1459(5) 


0.0848(5) 


0.0226(5) 


0.2087(4) 


0.1460(5) 


0.0848(5) 


0.0226(5) 


(75 i 




0.2203(6) 


0.1567(6) 


0.0948(5) 


0.0318(4) 


0.2203(5) 


0.1567(5) 


0.0947(5) 


0.0317(3) 


(75 1 


5£k) 


0.2268(6) 


0.1633(7) 


0.1021(7) 


0.0393(4) 


0.2268(6) 


0.1634(6) 


0.1021(7) 


0.0392(4) 


(15 i 


5/) 


0.2324(8) 


0.170(1) 


0.110(1) 


0.048(1) 


0.2325(7) 


0.1699(9) 


0.110(1) 


0.048(1) 




(b) 13 = 6.2 








Gauge invariant 






Non-invariant 




Operator 


m q a = 0.023 


m q a = 0.015 


m q a = 0.008 


m q a —> 


m q a = 0.023 


m q a = 0.015 


m q a = 0.008 


m q a — > 


(75 S 




0.0927(3) 


0.0604(2) 


0.0326(3) 


0.0004(4) 










(75 i 




0.1017(3) 


0.0679(3) 


0.0394(3) 


0.0058(4) 


0.1017(3) 


0.0679(3) 


0.0393(3) 


0.0058(4) 


(75 t 




0.1046(4) 


0.0706(3) 


0.0420(4) 


0.0083(4) 


0.1046(4) 


0.0706(3) 


0.0420(4) 


0.0083(4) 


(75 t 


3 €4) 


0.1062(4) 


0.0724(3) 


0.0438(4) 


0.0102(4) 


0.1063(4) 


0.0723(3) 


0.0438(4) 


0.0102(4) 


(75 i 




0.1019(3) 


0.0680(3) 


0.0393(3) 


0.0057(4) 


0.1021(4) 


0.0681(3) 


0.0394(3) 


0.0056(3) 


(75 i 




0.1047(3) 


0.0706(3) 


0.0419(3) 


0.0081(4) 


0.1049(4) 


0.0706(3) 


0.0420(4) 


0.0081(4) 


(75 t 




0.1064(4) 


0.0724(3) 


0.0438(3) 


0.0101(4) 


0.1066(4) 


0.0724(3) 


0.0439(4) 


0.0100(4) 


(75 i 


3/) 


0.1080(4) 


0.0742(3) 


0.0461(4) 


0.0127(4) 


0.1082(4) 


0.0742(4) 


0.0461(5) 


0.0125(4) 
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TABLE VI. Bare pion decay constant f£ in lattice units. The bottom line shows results 
obtained from pion operator fi P ^ 5 . 

(a) (3 = 6.0 



Gauge invariant Non-invariant 



Operator 


m q a = 0.030 


m q a = 0.020 


m q a = 0.010 


m q a — > 


m q a = 0.030 


m q a = 0.020 


m q a = 0.010 


m q a — > 


(7475 


®€s) 


0.0770(5) 


0.0686(6) 


0.0586(4) 


0.0495(6) 


0.0859(6) 


0.0772(6) 


0.0673(6) 


0.0582(8) 


(7475 


®£fc£s) 


0.0606(8) 


0.0551(6) 


0.0482(4) 


0.0420(5) 


0.081(1) 


0.0743(9) 


0.0659(6) 


0.0581(7) 


(7475 


® €fe?4) 


0.0481(8) 


0.0440(7) 


0.0389(4) 


0.0342(5) 


0.078(1) 


0.071(1) 


0.0640(7) 


0.0572(9) 


(7475 


®&) 


0.0390(8) 


0.0357(8) 


0.0308(4) 


0.0267(6) 


0.076(2) 


0.069(2) 


0.0613(7) 


0.0541(9) 


(7475 


®€4?5) 


0.093(2) 


0.084(2) 


0.073(1) 


0.064(1) 


0.093(2) 


0.084(2) 


0.073(1) 


0.064(1) 


(7475 


® &£m) 


0.076(1) 


0.069(1) 


0.0603(8) 


0.0524(9) 


0.086(1) 


0.079(1) 


0.0684(8) 


0.060(1) 


(7475 


®£fc) 


0.061(1) 


0.056(1) 


0.0477(6) 


0.0410(7) 


0.082(2) 


0.075(1) 


0.0646(8) 


0.0559(9) 


(7475 


®/) 


0.049(1) 


0.045(1) 


0.038(1) 


0.033(1) 


0.080(2) 


0.072(2) 


0.062(2) 


0.054(2) 


(75 ® 


&) 


0.0789(6) 


0.0697(6) 


0.0602(5) 


0.0509(6) 












(b) (3 = 6.2 








Gauge invariant 






Non-invariant 




Operator 


m q a = 0.023 


m q a = 0.015 


m q a = 0.008 


m q a — > 


m q a = 0.023 


m q a = 0.015 


m q a = 0.008 


m q a — » 


(7475 


®&) 


0.0524(8) 


0.0454(4) 


0.0404(4) 


0.0341(6) 


0.0594(6) 


0.0520(5) 


0.0468(6) 


0.040(1) 


(7475 ® 


0.044(1) 


0.0381(6) 


0.0344(5) 


0.0294(6) 


0.058(1) 


0.0511(7) 


0.0463(6) 


0.040(1) 


(7475 ® £kM 


0.0362(9) 


0.0315(6) 


0.0284(4) 


0.0243(5) 


0.058(1) 


0.0502(8) 


0.0455(6) 


0.0391(9) 


(7475 


®£ 4 ) 


0.031(1) 


0.0263(5) 


0.0235(5) 


0.0195(5) 


0.059(2) 


0.050(1) 


0.045(1) 


0.0376(7) 


(7475 


®«46) 


0.068(3) 


0.057(1) 


0.0509(6) 


0.043(1) 


0.068(3) 


0.057(1) 


0.0508(6) 


0.043(1) 


(7475 


® &£m) 


0.057(2) 


0.048(1) 


0.0425(6) 


0.035(1) 


0.064(3) 


0.054(1) 


0.0479(7) 


0.040(1) 


(7475 


®£fc) 


0.047(2) 


0.0390(9) 


0.0346(5) 


0.029(1) 


0.062(3) 


0.052(1) 


0.0461(7) 


0.038(2) 


(7475 


® J) 


0.039(2) 


0.0322(7) 


0.0285(5) 


0.024(1) 


0.062(3) 


0.051(1) 


0.0452(7) 


0.037(2) 


(75 ® 


&) 


0.055(1) 


0.0485(4) 


0.0425(6) 


0.036(1) 











TABLE VII. Perturbatively renormalized pion decay constants Z\ f£ in lattice unit. 



(a) (3 = 6.0 

Gauge invariant Non-invariant 



Operator 


m q a = 0.030 


m q a = 0.020 


m q a = 0.010 


m q a — > 


m q a = 0.030 


m q a = 0.020 


m q a = 0.010 


m q a — » 


(7475 ® £5) 


0.0770(5) 


0.0686(6) 


0.0586(4) 


0.0495(6) 


0.0766(5) 


0.0688(6) 


0.0600(5) 


0.0519(7) 


(7475 ® kk&) 


0.0693(9) 


0.0630(7) 


0.0551(5) 


0.0480(6) 


0.0695(9) 


0.0635(7) 


0.0563(5) 


0.0497(6) 


(7475 ® Zk£.i) 


0.066(1) 


0.0605(9) 


0.0534(5) 


0.0471(7) 


0.066(1) 


0.0609(9) 


0.0548(6) 


0.0490(8) 


(7475 ® £,4) 


0.058(1) 


0.053(1) 


0.0461(6) 


0.0399(9) 


0.065(1) 


0.059(1) 


0.0525(6) 


0.0463(8) 


(7475 ® ?4?5) 


0.074(2) 


0.067(1) 


0.0580(8) 


0.0502(8) 


0.073(1) 


0.067(1) 


0.0580(8) 


0.0503(9) 


(7475 ® klkm) 


0.071(1) 


0.065(1) 


0.0560(7) 


0.0487(8) 


0.071(1) 


0.065(1) 


0.0566(7) 


0.0493(8) 


(7475 ® kk) 


0.070(2) 


0.064(1) 


0.0546(7) 


0.0469(8) 


0.070(1) 


0.064(1) 


0.0552(7) 


0.0478(8) 


(7475 ® I) 


0.068(2) 


0.062(1) 


0.052(2) 


0.045(2) 


0.068(2) 


0.062(1) 


0.053(2) 


0.047(2) 




(b) 13 = 6.2 






Gauge invariant 






Non-invariant 




Operator 


m q a = 0.023 


m q a = 0.015 


m q a = 0.008 


m q a 


m q a = 0.023 


m q a = 0.015 


m q a = 0.008 


m q a — * 


(7475 ® 6) 


0.0524(8) 


0.0454(4) 


0.0404(4) 


0.0341(6) 


0.0530(6) 


0.0464(4) 


0.0417(5) 


0.0355(9) 


(7475 ® ffcfo) 


0.049(1) 


0.0432(7) 


0.0390(6) 


0.0334(7) 


0.051(1) 


0.0442(6) 


0.0400(5) 


0.0344(9) 


(7475 ® ik£.4.) 


0.049(1) 


0.0424(7) 


0.0381(5) 


0.0326(7) 


0.050(1) 


0.0435(7) 


0.0393(5) 


0.0338(8) 


(7475 ® 54) 


0.045(2) 


0.0383(8) 


0.0342(8) 


0.0284(8) 


0.051(2) 


0.0433(9) 


0.0390(8) 


0.0326(6) 


(7475 ® ?4&) 


0.055(2) 


0.046(1) 


0.0411(5) 


0.0344(9) 


0.055(2) 


0.046(1) 


0.0410(5) 


0.034(1) 


(7475 ® itkm) 


0.054(2) 


0.045(1) 


0.0398(6) 


0.033(1) 


0.054(2) 


0.0453(9) 


0.0402(5) 


0.034(1) 


(7475 ®ik) 


0.053(2) 


0.044(1) 


0.0392(6) 


0.033(2) 


0.054(2) 


0.0450(1) 


0.0399(6) 


0.033(2) 


(7475 ® I) 


0.052(3) 


0.043(1) 


0.0385(7) 


0.032(2) 


0.054(3) 


0.0444(9) 


0.0393(6) 


0.033(2) 
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TABLE VIII. Nonerturbatively renormalized pion decay constants Z A f£ in lattice unit. 











(a) 


1 /3 = 6.0 
















Gauge invariant 






Non-invariant 




Operator 


m g a = 0.030 


m q a = 0.020 


m q a = 0.010 


m q a — + 


m q a = 0.030 


m q a = 0.020 


m q a = 0.010 


m q a — » 


(7475 t 


Kb) 


0.0770(5) 


0.0686(6) 


0.0586(4) 


0.0495(6) 


0.0731(5) 


0.0656(5) 


0.0572(5) 


0.0494(7) 


(7475 t 




0.0728(9) 


0.0661(8) 


0.0579(5) 


0.0504(6) 


0.0693(9) 


0.0633(8) 


0.0562(5) 


0.0495(6) 


(7475 15 




0.071(1) 


0.065(1) 


0.0575(6) 


0.0507(8) 


0.067(1) 


0.0616(9) 


0.0554(6) 


0.0495(8) 


(7475 t 


554) 


0.071(1) 


0.065(1) 


0.0562(8) 


0.049(1) 


0.066(1) 


0.060(1) 


0.0535(6) 


0.0472(8) 


(7475 t 




0.074(2) 


0.067(1) 


0.0585(8) 


0.0507(8) 


0.074(2) 


0.067(1) 


0.0585(8) 


0.0507(9) 


(7475 t 




0.075(1) 


0.068(1) 


0.0594(8) 


0.0516(9) 


0.073(1) 


0.067(1) 


0.0582(7) 


0.0507(8) 


(7475 t 




0.075(2) 


0.068(1) 


0.0586(8) 


0.0505(8) 


0.072(1) 


0.066(1) 


0.0566(7) 


0.0490(8) 


(7475 t 




0.074(2) 


0.068(2) 


0.057(2) 


0.050(2) 


0.070(2) 


0.064(1) 


0.055(2) 


0.048(2) 




(b) 13 = 6.2 








Gauge invariant 






Non-invariant 




Operator 


m q a = 0.023 


m q a = 0.015 


m q a = 0.008 


m q a —* 


m q a = 0.023 


m q a = 0.015 


m q a = 0.008 


m q a — + 


(7475 C 




0.0524(8) 


0.0454(4) 


0.0404(4) 


0.0341(6) 


0.0513(6) 


0.0450(4) 


0.0404(5) 


0.0344(9) 


(7475 t 




0.051(1) 


0.0449(8) 


0.0405(6) 


0.0347(7) 


0.051(1) 


0.0441(6) 


0.0400(5) 


0.0344(9) 


(7475 


5 &&) 


0.052(1) 


0.0448(8) 


0.0403(5) 


0.0345(7) 


0.050(1) 


0.0439(7) 


0.0397(5) 


0.0342(8) 


(7475 G 


5?4) 


0.053(2) 


0.0451(9) 


0.0403(9) 


0.0335(9) 


0.052(2) 


0.0440(9) 


0.0396(8) 


0.0331(6) 


(7475 c 


5^5) 


0.055(2) 


0.047(1) 


0.0414(5) 


0.035(1) 


0.055(2) 


0.046(1) 


0.0413(5) 


0.035(1) 


(7475 c 


5 &£m) 


0.056(2) 


0.047(1) 


0.0418(6) 


0.035(1) 


0.055(2) 


0.046(1) 


0.0412(6) 


0.034(1) 


(7475 


5£k) 


0.056(3) 


0.047(1) 


0.0415(7) 


0.034(2) 


0.055(2) 


0.046(1) 


0.0407(6) 


0.034(2) 


(7475 


5^) 


0.056(3) 


0.047(1) 


0.0412(7) 


0.034(2) 


0.055(3) 


0.045(1) 


0.0401(7) 


0.033(2) 



TABLE IX. Pion mass squared (m^) 2 in GeV 2 . 







Gauge invariant 






Non-invariant 




Operator 


= 6.0 


= 6.2 


a -> 


= 6.0 


= 6.2 


a -> 


(75 ® 6) 


0.0066(9) 


0.003(3) 


0.000(6) 








(75 ® &&) 


0.085(2) 


0.042(4) 


0.002(8) 


0.085(2) 


0.042(4) 


0.002(8) 


(75 ® ZkU) 


0.118(3) 


0.061(4) 


0.006(9) 


0.118(3) 


0.061(4) 


0.005(9) 


(75 ® U) 


0.147(6) 


0.074(5) 


0.000(10) 


0.147(6) 


0.074(5) 


0.000(10) 


(75 ® £460 


0.084(2) 


0.042(3) 


0.001(6) 


0.084(2) 


0.041(3) 


0.000(6) 


(75 ® CeCm) 


0.118(3) 


0.059(4) 


0.002(8) 


0.118(3) 


0.059(4) 


0.002(8) 


(75 ® 60 


0.146(4) 


0.074(4) 


0.004(9) 


0.146(4) 


0.073(5) 


0.000(10) 


(75 ® /) 


0.178(5) 


0.092(4) 


0.009(9) 


0.178(6) 


0.091(5) 


0.010(10) 
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TABLE X. Perturbatively renormalized pion decay constants Z A in MeV. The bottom 

(P)5 

line shows results obtained with the pion operator in the £5 channel 







Gauge invariant 






Non- invariant 




Operator 


P — D.U 


a co 

p — o.z 


a — > u 


p — D.U 


a an 

p — o.z 


— * u 


(7475 ® £5) 


t)0[Z ) 


no/ q\ 
i)Z[6 ) 


QO/£\ 

oy(o J 


lUU^z J 


yo(4 j 


OQ/7\ 

92(7) 


(7475 ® ?fc?5j 


no/ 1 \ 

yz^l J 


on/ q\ 
i)U[o ) 


oo(0 J 


ofi/o\ 
»D(z J 


09/ 0\ 

i)o[o ) 


VU{D ) 


(7475 ® 6^4) 


91(2) 


88(3) 


85(5) 


94(2) 


91(2) 


88(5) 


(7475 O £4) 


77(2) 


77(3) 


77(5) 


89(2) 


88(2) 


87(4) 


(7475 ® 


97(2) 


93(3) 


89(7) 


97(2) 


93(4) 


88(7) 


(7475 O 


94(2) 


90(3) 


86(7) 


95(2) 


91(4) 


86(7) 


(7475 ®4) 


90(2) 


88(4) 


85(9) 


92(2) 


89(5) 


87(9) 


(7475 O J) 


87(4) 


86(5) 


85(10) 


90(4) 


88(5) 


86(11) 


(75 ® 6) 


98(1) 


94(3) 


89(6) 









TABLE XI. Non-perturbatively renormalized pion decay constants Z A in MeV unit. The 
bottom line for fn is reproduced from Table |X| for convenience. 







Gauge invariant 






Non-invariant 




Operator 


/3 = 6.0 


= 6.2 


a -> 


(3 = 6.0 


/3 = 6.2 


a -» 


(7475 ® £5) 


95(2) 


92(3) 


89(6) 


95(2) 


93(3) 


91(7) 


(7475 ® Cfc^s) 


97(2) 


94(3) 


90(6) 


95(2) 


93(3) 


90(6) 


(7475 ® 


98(2) 


93(3) 


89(6) 


95(2) 


92(2) 


89(5) 


(7475 ® £4) 


94(2) 


90(3) 


87(6) 


91(2) 


89(2) 


88(4) 


(7475 8> 


98(2) 


94(3) 


90(7) 


98(2) 


93(4) 


89(7) 


(7475 &> &£m) 


99(2) 


94(3) 


89(7) 


98(2) 


93(4) 


88(8) 


(7475 8> £fc) 


97(2) 


93(5) 


89(9) 


94(2) 


91(5) 


88(9) 


(7475 ® -0 


96(5) 


92(5) 


89(11) 


92(4) 


90(5) 


88(11) 


(75 ® 6) 


98(1) 


94(3) 


89(6) 
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